reset

a_96  =  2.956703416589903e+01;
a_97  =  3.538176138679367e+01;
a_98  =  4.519514157602768e+01;
a_99  =  6.760661759397387e+01;


l_96  = a_96**(-2./3); 
l_97  = a_97**(-2./3); 
l_98  = a_98**(-2./3); 
l_99  = a_99**(-2./3); 

set title 'Comparison of 3d simulation and radially symmetric case\n for cubic NLS'
set term postscript eps enhanced color
set output 'proof_of_convergence.pdf'

set xlabel 'r, radius'
set ylabel '|Psi|'
plot[0:0.4] 'check_0050.bin' t 'rad symm' w l, \
	    'g_0075.dat' t '3d simulation' w lp


set output
set term x11
